Genomic and phenotypic insights from an atlas of genetic effects on DNA methylation.

Characterizing genetic influences on DNA methylation (DNAm) provides an opportunity to understand mechanisms underpinning gene regulation and disease. In the present study, we describe results of DNAm quantitative trait locus (mQTL) analyses on 32,851 participants, identifying genetic variants associated with DNAm at 420,509 DNAm sites in blood. We present a database of >270,000 independent mQTLs, of which 8.5% comprise long-range (trans) associations. Identified mQTL associations explain 15-17% of the additive genetic variance of DNAm. We show that the genetic architecture of DNAm levels is highly polygenic. Using shared genetic control between distal DNAm sites, we constructed networks, identifying 405 discrete genomic communities enriched for genomic annotations and complex traits. Shared genetic variants are associated with both DNAm levels and complex diseases, but only in a minority of cases do these associations reflect causal relationships from DNAm to trait or vice versa, indicating a more complex genotype-phenotype map than previously anticipated.

Identification of loci where DNA methylation potentially mediates genetic risk of type 1 diabetes.

The risk of Type 1 Diabetes (T1D) comprises both genetic and environmental components. We investigated whether genetic susceptibility to T1D could be mediated by changes in DNA methylation, an epigenetic mechanism that potentially plays a role in autoimmune diabetes. From enrichment analysis, we found that there was a common genetic influence for both DNA methylation and T1D across the genome, implying that methylation could be either on the causal pathway to T1D or a non-causal biomarker of T1D genetic risk. Using data from a general population comprising blood samples taken at birth (n = 844), childhood (n = 846) and adolescence (n = 907), we then evaluated the associations between 64 top GWAS single nucleotide polymorphisms (SNPs) and DNA methylation levels at 55 non-HLA loci. We identified 95 proximal SNP-cytosine phosphate guanine (CpG) pairs (cis) and 1 distal SNP-CpG association (trans) consistently at birth, childhood, and adolescence. Combining genetic co-localization and Mendelian Randomization analysis, we provided evidence that at 5 loci, ITGB3BP, AFF3, PTPN2, CTSH and CTLA4, DNA methylation is potentially mediating the genetic risk of T1D mainly by influencing local gene expression.

The epigenetic clock and physical development during childhood and adolescence: longitudinal analysis from a UK birth cohort.

Background: Statistical models that use an individual's DNA methylation levels to estimate their age (known as epigenetic clocks) have recently been developed, with 96% correlation found between epigenetic and chronological age. We postulate that differences between estimated and actual age [age acceleration (AA)] can be used as a measure of developmental age in early life. Methods: We obtained DNA methylation measures at three time points (birth, age 7 years and age 17 years) in 1018 children from the Avon Longitudinal Study of Parents and Children (ALSPAC). Using an online calculator, we estimated epigenetic age, and thus AA, for each child at each time point. We then investigated whether AA was prospectively associated with repeated measures of height, weight, body mass index (BMI), bone mineral density, bone mass, fat mass, lean mass and Tanner stage. Results: Positive AA at birth was associated with higher average fat mass [1321 g per year of AA, 95% confidence interval (CI) 386, 2256 g] from birth to adolescence (i.e. from age 0-17 years) and AA at age 7 was associated with higher average height (0.23 cm per year of AA, 95% CI 0.04, 0.41 cm). Conflicting evidence for the role of AA (at birth and in childhood) on changes during development was also found, with higher AA being positively associated with changes in weight, BMI and Tanner stage, but negatively with changes in height and fat mass. Conclusions: We found evidence that being ahead of one's epigenetic age acceleration is related to developmental characteristics during childhood and adolescence. This demonstrates the potential for using AA as a measure of development in future research.

Collapsed methylation quantitative trait loci analysis for low frequency and rare variants.

BACKGROUND: Single variant approaches have been successful in identifying DNA methylation quantitative trait loci (mQTL), although as with complex traits they lack the statistical power to identify the effects from rare genetic variants. We have undertaken extensive analyses to identify regions of low frequency and rare variants that are associated with DNA methylation levels. METHODS: We used repeated measurements of DNA methylation from five different life stages in human blood, taken from the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort. Variants were collapsed across CpG islands and their flanking regions to identify variants collectively associated with methylation, where no single variant was individually responsible for the observed signal. All analyses were undertaken using the sequence kernel association test. RESULTS: For loci where no individual variant mQTL was observed based on a single variant analysis, we identified 95 unique regions where the combined effect of low frequency variants (MAF ≤ 5%) provided strong evidence of association with methylation. For loci where there was previous evidence of an individual variant mQTL, a further 3 regions provided evidence of association between multiple low frequency variants and methylation levels. Effects were observed consistently across 5 different time points in the lifecourse and evidence of replication in the TwinsUK and Exeter cohorts was also identified. CONCLUSION: We have demonstrated the potential of this novel approach to mQTL analysis by analysing the combined effect of multiple low frequency or rare variants. Future studies should benefit from applying this approach as a complementary follow up to single variant analyses.

Systematic identification of genetic influences on methylation across the human life course.

BACKGROUND: The influence of genetic variation on complex diseases is potentially mediated through a range of highly dynamic epigenetic processes exhibiting temporal variation during development and later life. Here we present a catalogue of the genetic influences on DNA methylation (methylation quantitative trait loci (mQTL)) at five different life stages in human blood: children at birth, childhood, adolescence and their mothers during pregnancy and middle age. RESULTS: We show that genetic effects on methylation are highly stable across the life course and that developmental change in the genetic contribution to variation in methylation occurs primarily through increases in environmental or stochastic effects. Though we map a large proportion of the cis-acting genetic variation, a much larger component of genetic effects influencing methylation are acting in trans. However, only 7 % of discovered mQTL are trans-effects, suggesting that the trans component is highly polygenic. Finally, we estimate the contribution of mQTL to variation in complex traits and infer that methylation may have a causal role consistent with an infinitesimal model in which many methylation sites each have a small influence, amounting to a large overall contribution. CONCLUSIONS: DNA methylation contains a significant heritable component that remains consistent across the lifespan. Our results suggest that the genetic component of methylation may have a causal role in complex traits. The database of mQTL presented here provide a rich resource for those interested in investigating the role of methylation in disease.

DNA Methylation and BMI: Investigating Identified Methylation Sites at HIF3A in a Causal Framework.

Multiple differentially methylated sites and regions associated with adiposity have now been identified in large-scale cross-sectional studies. We tested for replication of associations between previously identified CpG sites at HIF3A and adiposity in ∼1,000 mother-offspring pairs from the Avon Longitudinal Study of Parents and Children (ALSPAC). Availability of methylation and adiposity measures at multiple time points, as well as genetic data, allowed us to assess the temporal associations between adiposity and methylation and to make inferences regarding causality and directionality. Overall, our results were discordant with those expected if HIF3A methylation has a causal effect on BMI and provided more evidence for causality in the reverse direction (i.e., an effect of BMI on HIF3A methylation). These results are based on robust evidence from longitudinal analyses and were also partially supported by Mendelian randomization analysis, although this latter analysis was underpowered to detect a causal effect of BMI on HIF3A methylation. Our results also highlight an apparent long-lasting intergenerational influence of maternal BMI on offspring methylation at this locus, which may confound associations between own adiposity and HIF3A methylation. Further work is required to replicate and uncover the mechanisms underlying the direct and intergenerational effect of adiposity on DNA methylation.

Prenatal and early life influences on epigenetic age in children: a study of mother-offspring pairs from two cohort studies.

DNA methylation-based biomarkers of aging are highly correlated with actual age. Departures of methylation-estimated age from actual age can be used to define epigenetic measures of child development or age acceleration (AA) in adults. Very little is known about genetic or environmental determinants of these epigenetic measures of aging. We obtained DNA methylation profiles using Infinium HumanMethylation450 BeadChips across five time-points in 1018 mother-child pairs from the Avon Longitudinal Study of Parents and Children. Using the Horvath age estimation method, we calculated epigenetic age for these samples. AA was defined as the residuals from regressing epigenetic age on actual age. AA was tested for associations with cross-sectional clinical variables in children. We identified associations between AA and sex, birth weight, birth by caesarean section and several maternal characteristics in pregnancy, namely smoking, weight, BMI, selenium and cholesterol level. Offspring of non-drinkers had higher AA on average but this difference appeared to resolve during childhood. The associations between sex, birth weight and AA found in ARIES were replicated in an independent cohort (GOYA). In children, epigenetic AA measures are associated with several clinically relevant variables, and early life exposures appear to be associated with changes in AA during adolescence. Further research into epigenetic aging, including the use of causal inference methods, is required to better our understanding of aging.

Genome-wide association studies identify genetic loci for low von Willebrand factor levels.

Low von Willebrand factor (VWF) levels are associated with bleeding symptoms and are a diagnostic criterion for von Willebrand disease, the most common inherited bleeding disorder. To date, it is unclear which genetic loci are associated with reduced VWF levels. Therefore, we conducted a meta-analysis of genome-wide association studies to identify genetic loci associated with low VWF levels. For this meta-analysis, we included 31 149 participants of European ancestry from 11 community-based studies. From all participants, VWF antigen (VWF:Ag) measurements and genome-wide single-nucleotide polymorphism (SNP) scans were available. Each study conducted analyses using logistic regression of SNPs on dichotomized VWF:Ag measures (lowest 5% for blood group O and non-O) with an additive genetic model adjusted for age and sex. An inverse-variance weighted meta-analysis was performed for VWF:Ag levels. A total of 97 SNPs exceeded the genome-wide significance threshold of 5 × 10(-8) and comprised five loci on four different chromosomes: 6q24 (smallest P-value 5.8 × 10(-10)), 9q34 (2.4 × 10(-64)), 12p13 (5.3 × 10(-22)), 12q23 (1.2 × 10(-8)) and 13q13 (2.6 × 10(-8)). All loci were within or close to genes, including STXBP5 (Syntaxin Binding Protein 5) (6q24), STAB5 (stabilin-5) (12q23), ABO (9q34), VWF (12p13) and UFM1 (ubiquitin-fold modifier 1) (13q13). Of these, UFM1 has not been previously associated with VWF:Ag levels. Four genes that were previously associated with VWF levels (VWF, ABO, STXBP5 and STAB2) were also associated with low VWF levels, and, in addition, we identified a new gene, UFM1, that is associated with low VWF levels. These findings point to novel mechanisms for the occurrence of low VWF levels.

Sixteen new lung function signals identified through 1000 Genomes Project reference panel imputation.

Lung function measures are used in the diagnosis of chronic obstructive pulmonary disease. In 38,199 European ancestry individuals, we studied genome-wide association of forced expiratory volume in 1 s (FEV1), forced vital capacity (FVC) and FEV1/FVC with 1000 Genomes Project (phase 1)-imputed genotypes and followed up top associations in 54,550 Europeans. We identify 14 novel loci (P<5 × 10(-8)) in or near ENSA, RNU5F-1, KCNS3, AK097794, ASTN2, LHX3, CCDC91, TBX3, TRIP11, RIN3, TEKT5, LTBP4, MN1 and AP1S2, and two novel signals at known loci NPNT and GPR126, providing a basis for new understanding of the genetic determinants of these traits and pulmonary diseases in which they are altered.

Multi-ancestry genome-wide association study of 21,000 cases and 95,000 controls identifies new risk loci for atopic dermatitis.

Genetic association studies have identified 21 loci associated with atopic dermatitis risk predominantly in populations of European ancestry. To identify further susceptibility loci for this common, complex skin disease, we performed a meta-analysis of >15 million genetic variants in 21,399 cases and 95,464 controls from populations of European, African, Japanese and Latino ancestry, followed by replication in 32,059 cases and 228,628 controls from 18 studies. We identified ten new risk loci, bringing the total number of known atopic dermatitis risk loci to 31 (with new secondary signals at four of these loci). Notably, the new loci include candidate genes with roles in the regulation of innate host defenses and T cell function, underscoring the important contribution of (auto)immune mechanisms to atopic dermatitis pathogenesis.

Lymphoblastoid cell lines reveal associations of adult DNA methylation with childhood and current adversity that are distinct from whole blood associations.

BACKGROUND: Some cohort studies bank lymphoblastoid cell lines (LCLs) as a renewable source of participant DNA. However, although LCL DNA has proved valuable for genetic studies, its utility in epigenetic epidemiology research is unknown. METHODS: To assess whether LCL DNA can be used for life-course environmental epigenomics, we carried out a pilot methylomic study (using the Illumina Infinium Human Methylation 450 BeadChip) of nil-passage, Epstein-Barr virus (EBV)-transformed LCLs (n = 42) and 28 matched whole-blood (WB) samples. These were from adult male participants of the British 1958 birth cohort, selected for extremes of social economic position (SEP) in childhood and adulthood, with additional information available on childhood abuse and prenatal tobacco exposure. RESULTS: We identified a small number of weak associations between these exposures and methylation levels of both individual CpG sites and genomic regions in WB and LCLs. However, only one of the regional, and none of the individual CpG site associations were common to both sample types. The lack of overlap between the associations detected in LCL compared with those found in WB could either be due to the EBV-transformation process, or to the fact that, unlike WB, LCLs are essentially a single (CD19+) cell type. We provide evidence that the latter is the more potent explanation, by showing that CpG sites known to be differentially methylated between different types of blood cell have significantly lower correlations (R = 0.11) than average (R = 0.2) between WB and LCLs in our datasets, whereas sites known to be affected by EBV-transformation have significantly higher correlations (R = 0.3). CONCLUSIONS: This small pilot study suggests that the DNA methylation profile of LCLs is more closely related to that of B cells than WB and, additionally, that LCLs may nevertheless be useful for life-course environmental epigenomics.

Longitudinal analysis of DNA methylation associated with birth weight and gestational age.

Gestational age (GA) and birth weight have been implicated in the determination of long-term health. It has been hypothesized that changes in DNA methylation may mediate these long-term effects. We obtained DNA methylation profiles from cord blood and peripheral blood at ages 7 and 17 in the same children from the Avon Longitudinal Study of Parents and Children. Repeated-measures data were used to investigate changes in birth-related methylation during childhood and adolescence. Ten developmental phenotypes (e.g. height) were analysed to identify possible mediation of health effects by DNA methylation. In cord blood, methylation at 224 CpG sites was found to be associated with GA and 23 CpG sites with birth weight. Methylation changed in the majority of these sites over time, but neither birth characteristic was strongly associated with methylation at age 7 or 17 (using a conservative correction for multiple testing of P < 1.03 × 10(-7)), suggesting resolution of differential methylation by early childhood. Associations were observed between birth weight-associated CpG sites and phenotypic characteristics in childhood. One strong association involved birth weight, methylation of a CpG site proximal to the NFIX locus and bone mineral density at age 17. Analysis of serial methylation from birth to adolescence provided evidence for a lack of persistence of methylation differences beyond early childhood. Sites associated with birth weight were linked to developmental genes and have methylation levels which are associated with developmental phenotypes. Replication and interrogation of causal relationships are needed to substantiate whether methylation differences at birth influence the association between birth weight and development.

Prenatal exposure to maternal smoking and offspring DNA methylation across the lifecourse: findings from the Avon Longitudinal Study of Parents and Children (ALSPAC).

Maternal smoking during pregnancy has been found to influence newborn DNA methylation in genes involved in fundamental developmental processes. It is pertinent to understand the degree to which the offspring methylome is sensitive to the intensity and duration of prenatal smoking. An investigation of the persistence of offspring methylation associated with maternal smoking and the relative roles of the intrauterine and postnatal environment is also warranted. In the Avon Longitudinal Study of Parents and Children, we investigated associations between prenatal exposure to maternal smoking and offspring DNA methylation at multiple time points in approximately 800 mother-offspring pairs. In cord blood, methylation at 15 CpG sites in seven gene regions (AHRR, MYO1G, GFI1, CYP1A1, CNTNAP2, KLF13 and ATP9A) was associated with maternal smoking, and a dose-dependent response was observed in relation to smoking duration and intensity. Longitudinal analysis of blood DNA methylation in serial samples at birth, age 7 and 17 years demonstrated that some CpG sites showed reversibility of methylation (GFI1, KLF13 and ATP9A), whereas others showed persistently perturbed patterns (AHRR, MYO1G, CYP1A1 and CNTNAP2). Of those showing persistence, we explored the effect of postnatal smoke exposure and found that the major contribution to altered methylation was attributed to a critical window of in utero exposure. A comparison of paternal and maternal smoking and offspring methylation showed consistently stronger maternal associations, providing further evidence for causal intrauterine mechanisms. These findings emphasize the sensitivity of the methylome to maternal smoking during early development and the long-term impact of such exposure.

Large-scale genome-wide association studies and meta-analyses of longitudinal change in adult lung function.

BACKGROUND: Genome-wide association studies (GWAS) have identified numerous loci influencing cross-sectional lung function, but less is known about genes influencing longitudinal change in lung function. METHODS: We performed GWAS of the rate of change in forced expiratory volume in the first second (FEV1) in 14 longitudinal, population-based cohort studies comprising 27,249 adults of European ancestry using linear mixed effects model and combined cohort-specific results using fixed effect meta-analysis to identify novel genetic loci associated with longitudinal change in lung function. Gene expression analyses were subsequently performed for identified genetic loci. As a secondary aim, we estimated the mean rate of decline in FEV1 by smoking pattern, irrespective of genotypes, across these 14 studies using meta-analysis. RESULTS: The overall meta-analysis produced suggestive evidence for association at the novel IL16/STARD5/TMC3 locus on chromosome 15 (P  =  5.71 × 10(-7)). In addition, meta-analysis using the five cohorts with ≥3 FEV1 measurements per participant identified the novel ME3 locus on chromosome 11 (P  =  2.18 × 10(-8)) at genome-wide significance. Neither locus was associated with FEV1 decline in two additional cohort studies. We confirmed gene expression of IL16, STARD5, and ME3 in multiple lung tissues. Publicly available microarray data confirmed differential expression of all three genes in lung samples from COPD patients compared with controls. Irrespective of genotypes, the combined estimate for FEV1 decline was 26.9, 29.2 and 35.7 mL/year in never, former, and persistent smokers, respectively. CONCLUSIONS: In this large-scale GWAS, we identified two novel genetic loci in association with the rate of change in FEV1 that harbor candidate genes with biologically plausible functional links to lung function.

Variability in the common genetic architecture of social-communication spectrum phenotypes during childhood and adolescence.

BACKGROUND: Social-communication abilities are heritable traits, and their impairments overlap with the autism continuum. To characterise the genetic architecture of social-communication difficulties developmentally and identify genetic links with the autistic dimension, we conducted a genome-wide screen of social-communication problems at multiple time-points during childhood and adolescence. METHODS: Social-communication difficulties were ascertained at ages 8, 11, 14 and 17 years in a UK population-based birth cohort (Avon Longitudinal Study of Parents and Children; N ≤ 5,628) using mother-reported Social Communication Disorder Checklist scores. Genome-wide Complex Trait Analysis (GCTA) was conducted for all phenotypes. The time-points with the highest GCTA heritability were subsequently analysed for single SNP association genome-wide. Type I error in the presence of measurement relatedness and the likelihood of observing SNP signals near known autism susceptibility loci (co-location) were assessed via large-scale, genome-wide permutations. Association signals (P ≤ 10-5) were also followed up in Autism Genetic Resource Exchange pedigrees (N = 793) and the Autism Case Control cohort (Ncases/Ncontrols = 1,204/6,491). RESULTS: GCTA heritability was strongest in childhood (h2(8 years) = 0.24) and especially in later adolescence (h2(17 years) = 0.45), with a marked drop during early to middle adolescence (h2(11 years) = 0.16 and h2(14 years) = 0.08). Genome-wide screens at ages 8 and 17 years identified for the latter time-point evidence for association at 3p22.2 near SCN11A (rs4453791, P = 9.3 × 10-9; genome-wide empirical P = 0.011) and suggestive evidence at 20p12.3 at PLCB1 (rs3761168, P = 7.9 × 10-8; genome-wide empirical P = 0.085). None of these signals contributed to risk for autism. However, the co-location of population-based signals and autism susceptibility loci harbouring rare mutations, such as PLCB1, is unlikely to be due to chance (genome-wide empirical Pco-location = 0.007). CONCLUSIONS: Our findings suggest that measurable common genetic effects for social-communication difficulties vary developmentally and that these changes may affect detectable overlaps with the autism spectrum.

Differences in smoking associated DNA methylation patterns in South Asians and Europeans.

BACKGROUND: DNA methylation is strongly associated with smoking status at multiple sites across the genome. Studies have largely been restricted to European origin individuals yet the greatest increase in smoking is occurring in low income countries, such as the Indian subcontinent. We determined whether there are differences between South Asians and Europeans in smoking related loci, and if a smoking score, combining all smoking related DNA methylation scores, could differentiate smokers from non-smokers. RESULTS: Illumina HM450k BeadChip arrays were performed on 192 samples from the Southall And Brent REvisited (SABRE) cohort. Differential methylation in smokers was identified in 29 individual CpG sites at 18 unique loci. Interaction between smoking status and ethnic group was identified at the AHRR locus. Ethnic differences in DNA methylation were identified in non-smokers at two further loci, 6p21.33 and GNG12. With the exception of GFI1 and MYO1G these differences were largely unaffected by adjustment for cell composition. A smoking score based on methylation profile was constructed. Current smokers were identified with 100% sensitivity and 97% specificity in Europeans and with 80% sensitivity and 95% specificity in South Asians. CONCLUSIONS: Differences in ethnic groups were identified in both single CpG sites and combined smoking score. The smoking score is a valuable tool for identification of true current smoking behaviour. Explanations for ethnic differences in DNA methylation in association with smoking may provide valuable clues to disease pathways.

No evidence for genome-wide interactions on plasma fibrinogen by smoking, alcohol consumption and body mass index: results from meta-analyses of 80,607 subjects.

Plasma fibrinogen is an acute phase protein playing an important role in the blood coagulation cascade having strong associations with smoking, alcohol consumption and body mass index (BMI). Genome-wide association studies (GWAS) have identified a variety of gene regions associated with elevated plasma fibrinogen concentrations. However, little is yet known about how associations between environmental factors and fibrinogen might be modified by genetic variation. Therefore, we conducted large-scale meta-analyses of genome-wide interaction studies to identify possible interactions of genetic variants and smoking status, alcohol consumption or BMI on fibrinogen concentration. The present study included 80,607 subjects of European ancestry from 22 studies. Genome-wide interaction analyses were performed separately in each study for about 2.6 million single nucleotide polymorphisms (SNPs) across the 22 autosomal chromosomes. For each SNP and risk factor, we performed a linear regression under an additive genetic model including an interaction term between SNP and risk factor. Interaction estimates were meta-analysed using a fixed-effects model. No genome-wide significant interaction with smoking status, alcohol consumption or BMI was observed in the meta-analyses. The most suggestive interaction was found for smoking and rs10519203, located in the LOC123688 region on chromosome 15, with a p value of 6.2 × 10(-8). This large genome-wide interaction study including 80,607 participants found no strong evidence of interaction between genetic variants and smoking status, alcohol consumption or BMI on fibrinogen concentrations. Further studies are needed to yield deeper insight in the interplay between environmental factors and gene variants on the regulation of fibrinogen concentrations.

Common variation contributes to the genetic architecture of social communication traits.

BACKGROUND: Social communication difficulties represent an autistic trait that is highly heritable and persistent during the course of development. However, little is known about the underlying genetic architecture of this phenotype. METHODS: We performed a genome-wide association study on parent-reported social communication problems using items of the children's communication checklist (age 10 to 11 years) studying single and/or joint marker effects. Analyses were conducted in a large UK population-based birth cohort (Avon Longitudinal Study of Parents and their Children, ALSPAC, N = 5,584) and followed-up within a sample of children with comparable measures from Western Australia (RAINE, N = 1364). RESULTS: Two of our seven independent top signals (P-discovery <1.0E-05) were replicated (0.009

Meta-analysis of genome-wide association studies identifies ten loci influencing allergic sensitization.

Allergen-specific immunoglobulin E (present in allergic sensitization) has a central role in the pathogenesis of allergic disease. We performed the first large-scale genome-wide association study (GWAS) of allergic sensitization in 5,789 affected individuals and 10,056 controls and followed up the top SNP at each of 26 loci in 6,114 affected individuals and 9,920 controls. We increased the number of susceptibility loci with genome-wide significant association with allergic sensitization from three to ten, including SNPs in or near TLR6, C11orf30, STAT6, SLC25A46, HLA-DQB1, IL1RL1, LPP, MYC, IL2 and HLA-B. All the top SNPs were associated with allergic symptoms in an independent study. Risk-associated variants at these ten loci were estimated to account for at least 25% of allergic sensitization and allergic rhinitis. Understanding the molecular mechanisms underlying these associations may provide new insights into the etiology of allergic disease.